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Abstract 

One of the key ideas in describing the multiparticle production 
at high energies in the perturbative QCD is the angular ordering in 
successive soft-gluon emission. We analyse the angular distribution 
of particles in a jet, and investigate the corrections to the angular 
ordering. At a small angle from the jet direction, the angular or- 
dering is exact at the next-to-leading order. At a large angle, the 
angular ordering overestimates the angular particle density by a fac- 
tor (1 — 7 In cos 2 (#/2)), where 7 is the anomalous dimension of the 
multiplicity. The 0(7) correction restores the boost invariance of the 
cross section. 
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1 Introduction 



The modified leading-log approximation (MLLA) |4[] is the formalism to 
systematically resum the double- and single-logarithms occurring in analysing 
the multiparticle production at high energies in the perturbative QCD. After 
the resummation, the perturbation theory is reorganised in power series of 7, 
where 7 is the anomalous dimension of the multiplicity: 7 ~ J^ot s / (27r). The 
hypothesis of the local parton-hadron duality (LPHD)000, on the other 
hand, predicts that the partonic particle flow at high energies (calculated 
in the perturbation theory) is proportional to the hadronic flow (to be ob- 
served in the experiments). This framework (MLLA+LPHD) had remarkable 
successes particularly in predicting the shape of the particle spectrum in the 
e + e~-annihilation, in excellent agreement with the experimental dataJHH ]/]]. 

One of the key ideas underlying the MLLA is the angular ordering (AO) M 
in the successive soft-gluon emission. Because of the presence of the soft 
singularity in QCD, the majority of the particles originate from the softQ 
gluons: An energetic parton emits soft gluons, the soft gluons in turn emit 
softer gluons, and so on. The cascade of soft gluons thus gives rise to large 
multiplicities at high energies. A large-angle soft gluon, however, cannot re- 
solve the transverse polarisation in a jet at scales smaller than the inverse 
of the transverse momentum of the soft gluon. Its emission, therefore, does 
not take place independently from each individual parton in a narrow jet. 
A large- angle soft gluon is emitted from a narrow jet as though from a sin- 
gle energetic particle with the net charge of the jet ( "coherence" )@| 10]. It 
implies that only those gluons that are emitted at small angles participate 
the cascade of partons and contribute to the major part of the multiplicity 
(AO). The coherence is also the origin of the characteristic bell-shape of the 
particle spectrum observed in the e + e~-annihilation[^]. 

One of the important applications of the MLLA is the coherent shower 
Monte Carlo simulation of the multiparticle production|]rTJ||12[. From the 



theoretical points of view, the simulation is meant to give the total particle 
multiplicity and the particle spectrum correctly (at the next-to-leading order 
of MLLA). The coherent shower Monte Carlo programs seem to be able to 

1 By "soft" we mean at much lower energies than the jet energy, but still at much 
higher than the inverse of the hadron size so that the effective QCD coupling for their 
emission is small 
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reproduce universally well the events in the e + e~-annihilation. 

While the AO is a powerful picture which enables us to organise the 
parton cascade as a simple Markov process, it is as well important to closely 
examine its limitation. The picture of coherent soft-gluon emission may 
appear to be a rough approximation which holds only when a soft gluon is 
emitted from a narrow jet. When the jet is not necessarily very collimated, 
a more quantitative analysis is required. It has been known that the AO 
is correct at the next-to-leading order (0(a s )) for angle-averaged quantities 
(exact angular ordering). For such quantities, therefore, the correction due to 
the transverse polarisation of a jet ("dipole correction") appears at the next- 
to-next order (0(7 2 )). This type of correction was evaluated at the lowest 
order (two-gluon emission cross section from qq, 0(a 2 s )) by Dokshitzer et 
in the planar gauge. They made use of the result to evaluate the 
correction to the total multiplicity. 

In the quantities which depend on the angles, on the other hand, the 
dipole correction can appear at the next-to-leading order c. In this article, 
we shall closely examine the angular distribution of particles in a jet. At the 
leading order (the double-log approximation), the direction of the registered 
particle may be identified with the direction of the subjet of the primary 
gluon (the gluon directly emitted from the high-energy parton) to which it 
belongs. Owing to the factorisation of the soft-gluon emission amplitude, 
the angular density of particles is given at this order as the product of two 
amplitudesfl3|1: the probability of the emission of a primary gluon at the 
angle of the registered particle, and the particle multiplicity in the subjet 
with an opening angle of the size of the emission angle of the primary gluon. 
This approximation is justified under the strong AO. 

At higher orders, we should take into account the difference in direction 
between the primary gluon and the registered particle. The density of the 
secondary gluon, though sharply peaked in the direction of the primary gluon, 
spreads inside the angular cone of the size of the emission angle of the primary 
gluon. One might therefore expect that the sharp leading-order peak of the 
distribution in the direction of a jet may be somewhat "smeared" by the 
spread of the secondary gluon. We shall find, however, that such smearing 
of the peak does not take place. Namely, the AO is exact at small angles. 
As for the particle density far from the jet direction, on the other hand, the 
AO has to be modified at the next-to-leading order. We shall find, however, 
that the correction is due to the lack of the Lorentz covariance in the angular 
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variables. If we formulate the MLLA in terms of Lorentz invariants, as in 



15 ] [14 1 , we do not have the dipole correction. 

In the next section, we review the AO formalism. Though the formulae 
presented in this section are all well-known, we include them in order to make 
the text self-contained. We first derive the AO for the angle-averaged cross 
section for one-gluon emission (0(a s )). The result is generalised to all orders 
in /alpha s . The AO is exact at the next-to-leading order (0(a s )). The AO, 
however, is by no means obvious for the quantities which depend explicitly 
on the angle. 

In sect. 3, we discuss how much the exact soft-gluon emission amplitudes 
are modified by the approximation of the AO. We analyse in detail the Feyn- 
man diagrams for two-gluon emission (at the tree level), and examine the 
possible next-to-leading correction in the one-particle-inclusive (1PI) cross 
section. The result is generalised in sect. 4 beyond the fixed order perturba- 
tion: i.e. the large logarithms are resummed to all orders in a s (MLLA). 

We find that there is no correction to the AO at a small angle. This is 
owing to the symmetry of the cross section under the exchange between the 
primary and the secondary gluons: After integrating the two-gluon emission 
cross section over the momentum of the primary gluon, we are left with the 
same angular form for the secondary gluon as that for the primary gluonP| 

In sect. 5, we evaluate the correction at a large angle, where the dipole 
correction appears at the next-to-leading order. We find that the correction 
reduces the angular density by a factor (1 + 7 In cos 2 (#/2)). The correction 
restores the invariance of the cross section under the boost in the jet direction. 



2 The Angular Ordering 



In this section, we review the AO in the successively softer gluon emission § H ||| 10] 
First, we consider one soft-gluon emission from the gg-pair. A simple form of 
the AO, which is exact at the next-to- leading order, is obtained after dividing 
the factorised soft-gluon emission amplitude into two terms and performing 
azimuthal integration in each term. Then we analyse the soft-gluon emission 
from a gluon. We generalise the result to the successive soft-gluon emission 
at an arbitrary order in a s . 



2 This is the case only for a massless quark. When the quark jet is from a heavy quark, 
there is a substantial correction to the AO at the forward direction. See ref fllQl . 
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2.1 One-Gluon Emission 



The lowest order amplitude for one soft-gluon emission from gg-pair (Fig.l) 
is essentially the same as the well-known soft-photon emission amplitude in 
QED, with the QCD coupling a s in place of the fine structure constant a. 
In the soft limit, the gluon emission amplitude factors out from the lowest 
order gg-pair production cross section o"o(e + e~ — > qq): 

. , _ 4:7ra s C F d 3 k 2pi ■ p 2 

CpOi s dkdQ l — cos9 PlP2 

2n 2 k (1 - cos6> fcpi )(l - cos0 kp2 ) °"° 

where pi, p 2 and k are the momenta of the quark, the antiquark and the 
gluon respectively. 9 kpi (&kp 2 ) is the angle between the directions of motion 
of the gluon and the quark (antiquark), 9 PlP2 is the angle between the quark 
and antiquark, and Q is the solid angle of k. The quark colour charge Cp 
is defined by '^{j a T a ) a p = CV5 a/ 3, where r a is the quark colour matrix. 

C F = (iV c 2 - 1)/2JV C for SU(iV c ) (C F = 4/3 with N c = 3). 
The angular factor in (1), 

(112) = l-cos6 PlP2 

(1 - cos0 fepi )(l - cos6» fcp2 ) 

is positive definite. It is singular both as 9 kpi — > and as 9k P2 — > 0. We can 
write it as the sum of two terms each of which is singular in one direction 
only. 



(1 I 2) = - 1(1 I 2) + — ^ + {1^2} (3) 

1 ' J 2\ K 1 } l-cos6 kpi 1-cosVJ 

In the first parenthesis on the rhs of (3), the singularity in (1 | 2) as Q kp2 — > 
is cancelled by the third term. When we integrate it over the azimuthal angle 
(around the direction of pi), making use of the relation 

cos 9k P2 = cos 9k pi cos 9 P1P2 + sin 9 kpi sin 9 P1P2 cos , (4) 

we find that it vanishes outside a cone 9 kpi > 9 P1P2 : 

~ i n f d< ^ = i n ®(9piP2-9kpi] 

1 - cos 9 kp2 J 1 - cos 9 kpi 

(5) 
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where Q(x) is the step function: Q(x) = for x < and Q(x) = 1 for x > 0. 

This result is the simplest example of the exact AO. The AO is called 
"exact" in the following sense: With the angular restriction 8 kpi < 9 prp2) 
we can simply use 1/(1 — cos6 kpi ) for the emission amplitude. Although 
the exact amplitude before the azimuthal integration is a more complicated 
quantity, which is identical to 1/(1 — cos6 kpi ) only at a very small angle, the 
approximation is exact in evaluating an angle-averaged quantity. 

2.2 Soft-Gluon Emission from a Gluon 

Let us next suppose that an energetic quark emits a gluon at a momentum 
k\ at a small angle 9 klPl <C I. Two diagrams contribute to the soft gluon 
emission (at a momentum k 2 ) from them (Fig. 2), if the invariant mass k 2 ■ 
Pi is much smaller than ki • p\. The factorised amplitude in the eikonal 
approximation is 



When the emission angle of the soft gluon 8 k2P1 is much larger than 9 kipi , we 
have 6 k2kl ~ 9 k2Pl , and 



Making use of the commutation relation of the colour matrix (7), we obtain 



The colour structure on the rhs of (9) implies that the large-angle soft gluon 
emission factorises as though it were emitted before the quark emits a parallel 
gluon. In general, a large-angle soft gluon is emitted coherently from a narrow 
jet: the gluon cannot resolve the polarisation in the jet, and is emitted with 
its net charge. Consequently, the secondary soft-gluon emission from the 
gluons emitted inside a narrow jet can occur only at a small angle. 




K K/fa pi Pi/ Pi 



(8) 



k 2 -h k 2 (l - cos0 fe2fcl ) k 2 -pi k 2 (l - cos0 fe2Pl ) 
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Let us now analyse the coherence and the AO in the same process on the 
cross section level. The soft-gluon emission amplitude factorises from the 
e + e~ — > qqg cross section: 



da(e + e -> qqgg) ~ dcr(e + e ^ qqg)^^ 1 



n 2p! • p 2 



(fc 2 -pi)(k 2 -p 2 ) 



2 ( (fc 2 • ki)(k 2 -Pi) (k 2 • fa)(k 2 • p 2 ) (k 2 ■ Pi)(k 2 ■ p 2 ) 
d 3 k 2 



da(e e — > ggg) 



2tt 2 A; 2 



x 



C A 

C F a s (p 1 | p 2 ) + — a s {(fa | pi) + (fa | p 2 ) - (pi | p 2 )} 



where 



/, i \ l-cos6» felPl 
(fa I Pi) ^ 



(&i I P2) 



(1 - cos6» fc2fcl )(l - cos6» fc2Pl ) 

1 - cos 9 klP2 

(1 - cos6 k2kl )(l - cose k2P2 ) 



(P1IP2) = j. \ C °*° P1P2 R r - (11) 

(1 - cos&fc 2 pJ(l - cosO k2P2 ) 

The gluon colour charge C A is defined by £ = C A 5 CC ' (C A = 3 for 

a,b 

SU(3)). The first term, proportional to GV, represents the emission from the 
gg-pair. The contribution proportional to Ca consists of three terms : 

^1 {( kl \ Pl ) + (fa \ p 2 ) - ( Pl \ p 2 )} , (12) 

each of which has an angular form similar to the emission from the qq. In 
(12), the singularity in the direction of the quark (pi) is cancelled between 
the first and the third terms, while the singularity in the direction of the 
antiquark (p 2 ) is cancelled between the second and the third term. Therefore, 
the angular factor (12) has the singularity only in the direction of the gluon 
(fa). It thus corresponds to the emission from the fci-gluon (mainly directed 
near the direction of fa). 



,(10) 
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Let us now assume that the fci-gluon is emitted at a small angle from 
the quark (9k lPl <C PlP2 — @k lP2 )- As we have seen above, the soft-gluon 
amplitude from a gluon vanishes when the emission angle of the soft gluon 
from the fci-gluon 6 k , zkl is much larger than 9 kipi . We, therefore, are inter- 
ested in the case that the soft-gluon emission angle k2kl is comparable to or 
smaller than 6 kipi . It implies that the angle between the soft gluon and the 
antiquark 9 k2P2 is nearly equal to 6 I P1P2 (~ 9 klP2 ). Consequently, (12) can be 
approximated as 

(12) = | !>,) + - -„ -. ^—}. (13) 

2 { i - cos e k2kl i - cos e k2P1 j 

After the azimuthal angle averaging, we obtain 

(is) — - ° Aa ; &(e klP1 - e k2kl ) (14) 

1 — COS V kl k 2 

Thus the AO is exact at small angles. 

The result obtained above can be generalised to higher orders. Let us 
consider successive soft-gluon emission with strongly ordered soft momenta 
Pi ^ k\ ^> &2 3> k% ^> • • •. We construct Feynman diagrams by successive 
soft-gluon insertion into less soft lines. At each step, the factorised soft-gluon 
emission amplitude (at a momentum k i+ i) can be divided into two types of 
contribution. One of them is the soft-gluon emission from the colour configu- 
ration under which the previous i-th gluon (at a momentum fcj) was emitted. 
The first term (proportional to Cf) in the factorised amplitude of the ki- 
gluon emission in (10) is its simplest example. This type of contribution 
takes exactly the same form as the i-th gluon emission. 

On the other hand, the new colour configuration caused by the emission 
of the i-th gluon gives rise to a new contribution. Suppose that the emission 
angle of the i-th gluon from the (i — l)-th gluon is small. There are three 
types of contribution which are relevant to this process. One is from the 
interference term between the emission from the i-th gluon and the (i — l)-th 
gluon (for i — 1, it is the first term in (12)). The second type is from the 
interference between the i-the gluon and the sum of all the possible soft-gluon 
insertion to the final particles except the i-th and (i — l)-th gluons. For % — 1, 
it is the second term in (12). For i > 2, the soft-gluon insertion can be done 
on more than one lines (not just on the p2-line in the second term of (12)). 
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Note that the angular factor is insensitive to the direction of the momentum 
of the line on which the soft gluon is inserted, as is the case in the second 
term on the rhs of (13). With the help of the Ward-Takahashi identity 



we can easily prove that the sum of all the possible insertion gives the same 
angular factor as the insertion to a single (non-soft) parton which forms a 
colour singlet with the i-th and (i — l)-th gluons. 

The third type is from the interference between the (i — l)-th gluon and 
the sum of all the possible soft-gluon insertion to the final particles except 
the i-th and (i — l)-th gluons. The contribution from this type of inter- 
ference includes the independent emission from (i — l)-th gluon mentioned 
earlier. We, however, have an additional contribution, associated with the 
non-commutativity of the non-Abelian gauge coupling. The latter can be 
analysed in a similar way as the second type (for % = 1, it is the third term 
in (12)). 

We thus obtain an angular factor 

C A a s I (fci_ a | h) + z ~ a z \ 1 • (15) 

I 1 - COS e h+lki 1 - COS flfei+xfci-i J 

The azimuthal averaging of (15) leads to the exact AO form 

(15) — T-^f e («Wi - «Wxfc) ( 16 ) 

1 - COS0 fc ., lfci 



The picture obtained above can be summarised as follows pl[|: The initial 
gg-pair emits a number of gluons independently (the primary gluons). Each 
primary gluon emits softer secondary gluons independently of one another 
at angles smaller than the emission angle of the primary gluon. Each of the 
secondary gluons emits softer gluons independently of one another at angles 
smaller than the emission angle of the secondary gluon, and so on. 

In this article, we are interested in the 1PI cross section. If we integrate 
it over the registered momentum and divide it by the total cross section, we 
obtain the total particle multiplicity per event. In the soft-gluon emission de- 
scribed above, the parton cascade ramifies into many branches. Let us take 
the branch to which the registered particle belongs. The emission angles 
along this branch is in the decreasing order, 9 P1P2 > 9k ipi > 9k 2 k l > 9k 3 k 2 
There are other branches corresponding to the independent emission from 
the parent partons. The integration over the momentum of the unregistered 
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particles gives rise to the logarithmic divergence. Factorisation takes place 
for each soft-gluon emission amplitude (either virtual or real), and the can- 
cellation of the singularity between real emission of the unregistered particle 



and virtual emission takes place, owing to the KLN theorem|20|. There are 



no residual large logarithms after the cancellation. We therefore can neglect 
the independent emission in the evaluation of the 1PI cross section at the 
next-to-leading order. 

If we go back to the idea of the coherence in the soft-gluon emission 
mentioned earlier in this subsection, we obtain an alternative picture of the 
angle-ordered branching process: The largest-angle soft gluon from the quark 
(similarly from the antiquark) is first emitted at an angle smaller than the 
angle between the gg-pair and with the colour charge of a quark (the net 
colour of the jet), because of the coherence. Then from the second-largest- 
angle soft gluon is emitted from the quark in a similar way, and so on. 
These gluons are angle-ordered, but not energy-ordered. (As a matter of 
fact, the largest-angle gluon tends to be softer than others]!!), because of 
the transverse momentum cutoff in the emission.) 

On the other hand, the emitted softer gluons emit softer gluons. The 
emission angles of the latter gluons do not exceed the emission angle of the 
gluon from which they are emitted, and so on. Along this side of branching, 
gluons are both angle-ordered and energy-ordered. 

It is on the latter picture that the coherent shower Monte Carlo simu- 
lation of the QCD multiparticle production is based On the cross 



section level, the two pictures are indeed equivalent. The "independent" 
emission in the former picture can be angularly ordered, just with dropping 
the Bose-statistical factor l/n\, where n is the number of soft gluons emitted 
independently from the parent parton. 



3 Modification due to the Azimuthal Aver- 
aging 

As is described in the previous section, the AO is exact at the next-to-leading 
order for angle-averaged quantities, such as the total multiplicity and the 
particle energy spectrum. The AO, however, is by no means obvious for the 
angular distribution of particles. 
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In order to illustrate how the azimuth-averaged amplitude (AO with a 
step function) differs from the exact amplitude, let us examine the sim- 
plest case^TJ: One soft-gluon emission amplitude from an energetic quark- 
antiquark (qq) pair. 

In fact, the order among angles is not a Lorentz invariant quantity. There- 
fore, the distribution of the soft gluons generated with the AO depends on 
the Lorentz frame in which the angles are defined. Suppose, for example, 
that the initial qq pair is produced at the rectangle from a gauge-invariant 
vertex (e.g. the electromagnetic current). According to the AO, the soft- 
gluon radiation angle from the quark has to be less than ir/2: namely in 
the quark-side hemisphere divided by the plane perpendicular to the quark 
direction. Similarly, the soft-gluon radiation is in the antiquark-side hemi- 
sphere divided by the plane perpendicular to the antiquark direction. We, 
therefore, have four solid-angle regions divided by two planes (Fig. 3). In the 
region I in Fig. 3, we have gluon radiation only from the quark, while in the 
region II we have radiation only from the antiquark. In the region III, the 
gluon can come both from the quark and from the antiquark. Finally, in the 
region IV, we see no soft-gluon radiation. 

Now let us boost these angular regions to the centre-of-mass system of the 
qq (we neglect the recoil from the soft-gluon radiation). We then find an an- 
gular region where no gluons are emitted (around the direction of the boost). 
Such a region does not exist in the original soft-gluon emission amplitude (1), 
which is azimuthal-symmetric in the cm. frame. 

This shortcoming of the AO in the angular distribution may be avoided 
at 0(a s ) if we define it in the centre-of-mass frame of qq. In this case, the 
azimuthal integration is trivial. The division of the amplitude into two parts 
is in this frame 



1 — cos 2 9k Pl 1 — cos 9k Pl 1 — cos 9 kp2 

Indeed, there is no larger angle than tt, and the the AO does not restrict the 
emission angle. 

Whenever the AO is non-trivial, however, the azimuth-averaged expres- 
sion (with a step function) is not identical to the exact cross section. In the 
following, we examine if there is any next-to-leading order correction to the 
exact AO in the angular distribution of particles. The quantity is obtained 
by integrating the 1PI cross section over the registered energy. 
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Let us first analyse the two soft-gluon emission amplitude at the tree level 
in the cm. frame. The cross section is given by 



, n 2p! • p 2 



(h -pi){ki -p 2 ) 



C F a s 



(k 2 -pi)(k 2 -p 2 ) 



| C A a s \ 2fci • pi | 2ki-p 2 2p 1 -p 2 



2 \(k 2 ■ h)(k 2 -pi) (k 2 • h)(k 2 • p 2 ) (k 2 • pi)(k 2 ■ p 2 ) 
d 3 h d 3 k 2 

X (2tt) 2 A;i (2n) 2 k 2 [ j 

where we have assumed p\,p 2 ^> ki ^> k 2 (see (10)). 

Suppose that k 2 is the registered momentum. We integrate (18) over k\ to 
obtain the 1PI cross section. Note that (18) is symmetric under the exchange 
k\ <-> k 2 . Therefore, it can be written as 



dcitree = aoCpa s 



2pi • p 2 



(k 2 -pi)(k 2 -p 2 ) 



C F a s 



2pi • p 2 



-pi)(fci -p 2 ) 




(fci • fc 2 )(fci -Pi) (h ■ k^iykx ■ p 2 ) (h • pi)(ki • p 2 ) 

x J^J*h_ (19) 
{2nyk 2 {27i) 2 k 1 V } 

We may use (19) both for k\ > k 2 and for k 2 > ki. 

When integrated over ki, the first term in the curly bracket on the rhs of 
(19), proportional to the quark colour charge Cjr, gives a logarithmic diver- 
gence, which is cancelled by the virtual correction at the vertex of the gg-pair 
production. Similarly, the integration of the second term (proportional to the 
gluon colour charge Ca) over k\ with k\ < k 2 gives a logarithmic singularity, 
which is cancelled by the virtual correction at the /c 2 -gluon emission vertex. 
No large logarithms remain after the cancellation. Because we are interested 
in the terms with large double-logarithms, we consider only the second term 
with k\ > k 2 . 

If the registered angle 9 = 9 k2P1 is small, we can make use of the approx- 
imation: 

1 — cos 9 1 — cos 9k 2P2 



1 - cos6> fclfc2 )(l - cos6> fclPl ) (1 - cos6> felfc2 )(l - cos6» fclP2 ) 
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1 COS $ p\p2 



(1 - cos0 fclPl )(l - cos6» fclP2 ) 
1 — cos 9 



(1 - cos6> fclfc2 )(l - cos6> felPl ) + 1 - cos9 kl k 2 l-cos9 klPl ^ ^ 

— * 

Its azimuthal integration over k\ gives the step function Q(9—9 kl k 2 ), similarly 
to (14). Namely, the exact AO takes place also for the opposite energy or- 
dering (we call it the "reverse exact AO"), because of the symmetry between 
(18) and (19). We thus obtain 

do- C F a s C A u s f Pl dki f 1 d(cos6> fclfc2 ) 

<7 - 



d(cos9)dk 2 n k 2 n(l — cos 9) Jk 2 ki Jcose 1 — cos^^ 

C F a s C A a s gi r 1 d(cos 9 klk . 2 ) 
7r/c 2 7r(l — cos9) k 2 Jcose 1 — cos9 kl k 2 

The expression (20) has the form of the product of two amplitudes: the 
one-gluon emission amplitude (primary gluon) at the registered angle 9, and 
the soft-gluon emission amplitude (from the primary gluon) integrated over 
the phase space within the angle-ordered cone. The latter is identical to the 
gluon multiplicity (at 0(a s )) in the subjet of the opening angle identical to 
9. 

In the next section, where we evaluate the angular distribution to all or- 
ders in a s , we shall make systematical use of the identification of the emission 
angle of the primary gluon with the registered angle 9. The reverse exact AO 
above implies that at small angles there is no next-to-leading order correction 
to this approximation. 



4 Resummation of Double-Logarithms 

Because the integration over the momentum of each unregistered gluon gives 
rise to a double-logarithm, higher order contributions in the perturbative 
series in a s are not negligible in the high-energy limit compared with the 
contributions at the lowest orders. We, therefore, have to sum the large 
double-logarithms up to all orders in a s to redefine the perturbative expan- 
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sion (resummation). The resummation of the leading double logarithm is 
done by the double-log approximation^ (DLA)[Q[|TJ . 

When many gluons are emitted (with a high power in a s ), the matrix 
element shows, in general, complicated dependence on their angular con- 
figuration. The amplitude, however, takes particularly simple form for the 
configuration in which we are interested in DLA. The largest logarithm at 
each order in a s comes from the phase space configuration in which all the 
momenta are strongly ordered: 

Px > h > k 2 > h > • • • > q (22) 

where px is the energy of the initial quark, kx is that of the primary soft 
gluon, k 2 is that of softer secondary gluon, etc., and q is the energy of the 
registered particle. The emission angles are also strongly ordered. Because 
of the angular ordering (or coherence), we have 

9k lPl > 9 k , 2kl ^ 9k 3 k 2 > • • • (23) 

where 9 kipi is the emission angle of the primary soft gluon from the direction 
of the quark, 9 k2kl is that of the softer secondary gluon from the direction of 
the primary gluon, etc. If 9 is the angle between the registered particle and 
the energetic quark, we have 

9 ~ 9 klPl . (24) 
The ^-distribution is given at the leading order by 
d < n > C F a s [Pi dkx r d d9 1 d(pi 



5(0 - £4 



d# 7T 2 J h J 9 1 

S(0- Ox) + £ J] ° A ' f h <l/,J+J I'"' <U ^+' (lo '+ l 

n=l j=l 

2C F a s fPi dk 



it 2 J kj+x J 9 j+1 



i=0 



'J^L f 1 ^M g (k 2 9 2 ) . (9x = 9 klPl , 9 t = 9 kiki _ x for i > 2) 
tx9 J k 



3 This approximation, including its single-logarithmic next-to-leading order corrections, 
is given as a modified version of the more conventional leading-log approximation [|l7| ; 
hence the name "modified leading- log approximation" (MLLA). 
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Use is made of (24) to obtain the last expression, where M g (k 2 6 2 ) is the 
multiplicity in a gluon jet with the energy k and the opening angle 9: 

= i + e ft ^ f r f 11 ■ * = *■ 9 » = 9 > < 26 > 

The series summation can be represented as an integral equation 

Ms(m ^ff ffu.M<t) + l (27) 

At asymptotic high energies, the equation (27) can be approximated by an 
homogeneous equation: 



M^^-^fff'^mel). (28) 
The solution of the homogeneous equation (28) is given by 

M g (k 2 e 2 ) = c(k 2 e 2 y 

where the normalisation constant C is not determined by the homogeneous 
equation (28). 

The next-to- leading order corrections to (25) come from three cases Q [0 |L8 



(i) When the strong energy ordering (22) does not hold at one of the 
branching points in the successive gluon emission: 



kl 0(1) . (30) 



h-i 

For the non-soft emission, the integration kernel has to be modified in (25)- 
(28): 

/ ^— - — ► / dxi ( , (x\ = — ) (for the emission from a quark) 

Jo ki Jo \xi 4/ px 
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where Nf is the number of active quark flavours. 

(ii) The virtual corrections at the gluon-emission vertex run the effective 
QCD coupling. Replacing a s (k 2 J (k± = A^sin^) for the QCD coupling on 
the rhs of the integral equation (28) (together with the replacement (31)), 
we obtain its asymptotic solution: 



M g (k\) 



C 



71 



exp 



2 7 o V / ln(A;i/A| Ci? ) 



(32) 




12tt 



With the running coupling, the anomalous dimension 7 is replaced by 

7o 7i 



= d lnM g (kl) 



+ 



<ilnk 2 ± ^ln(k 2 JA 2 QCD ) Hk 2 JA 2 QCD ) 



(33) 



The logarithmic integration of M g gives a factor I/7, in stead of a logarithm. 

(iii) When one of the soft-gluon emission angle is of the same order of 
magnitude as the preceding one. If kiki _J6 ki _ lki _ 2 = 0(1) (i > 3) (but the 

— * — * 

strong AO holds between k\ and k 2 , & kipi 3> 0K 2 ki), then we can make use of 

— # — # 

the approximation (24), and the azimuthal averaging over ki (around 
leads to the exact AO discussed in sect. 2. 

The next-to-leading correction to the AO, therefore, can occur only for 
the configuration: 



~ 1 = 0(1), 9 kipi + 9 k2kl = 9 k2P1 ^ 



e 



kipi 



e . 



(34) 



In this case, we have to use a more exact form for the soft-gluon emission. 

At a small angle (fi C 1), however, we can make use of the reverse exact 
AO discussed in the last section. The angular particle density is given by 
multiplying the tree-level two-gluon emission cross section (divided by the 
total cross section) with the multiplicity in the fe-subjet with an opening 
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angle identical to the emission angle of the fc2-gluon 9k 2 kx- Following the 
steps which has led to (21), we obtain 

d < n > 2C F a s 2C A a s 1 m dJfei dk 2 f 9 d9 k 



d9 it it 9 J h 

2C F a s C A a s 1 1 dkt „ r 



el TTJ kj s( 2 W 

M 9 (A;^ 2 ) 



7r 2tt ~f 2 6 J k\ 

= 2 -^] M M0 2 ), (35) 

where use is made of the value of 7, (29), to obtain the last expression, which 
is identical to the leading order result (cf. the last expression in (25)). 

Finally, let us discuss the coherent Monte Carlo simulation which gener- 
ates particles with the AO0@. It appears that the particle flow generated 
by the simulation is in sharp contrast with the picture of the reverse exact 
AO, discussed in the previous section and used to prove (35) above. In 
the simulation, the fci-gluon emitted at a large angle can emit a softer k%- 
gluon at a similar large angle to form a small angle (i.e. 9 klPl + Ok 2 ki = 9, 
9k ipi > @k 2 ki 3> 9). Such a process would not give any contribution under 
the picture of the reverse exact AO in (21). 

In fact, when 9 kipi ^> 9, the partial amplitude (20) gives a positive 
contribution for 6k 2 ki < @k lPl , while the contribution becomes negative for 
&k 2 k 1 > 0k lPl - The contributions cancel out one another in average, which 
results in the reverse exact AO in (21). On the other hand, the Monte 
Carlo simulation assigns a positive probability to the former case, and zero 
to the latter. Thus the positive probability for the large angle process is 
not cancelled. This leads to the overestimation of the particle density at the 
next-to-leading order in the forward direction. 

The simulation, however, somewhat underestimates the contribution from 
the small angle process in which 9 > 9 k2kl > 9/2, because of the AO. In fact, 
after integrating over the direction of ki, the errors in the simulation cancel 
out. The exactness of the angular distribution at the next-to-leading order is 
required at small angles by the fact that the small angle region contributes 
to the total multiplicity at the leading order. Because the AO is exact in 
the total multiplicity, as was discussed in sect. 2, the particle density in the 
direction of a jet generated from a massless quark in the simulation with the 
exact AO cannot have next-to-leading order correction. 
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5 Corrections at Large Angles 



Let us next analyse the correction to the AO at a large angle. Because the 
large-angle particles contribute to the total multiplicity at the next-to-leading 
order, it is possible that the correction to the AO appears at the next order 
at a large angle . 

We start with the two-soft-gluon emission cross section at the tree level 
(19). We are interested in the term proportional to the gluon colour charge 
C A : 

2p x -p 2 d 3 k 2 d 3 h 

da(g^g) = a C F a s 



(k 2 ■ Pl )(k 2 ■ p 2 ) (2ir) 2 k 2 (2ir) 2 k 



i 



C A a s { 2k 2 ■ pi 2k 2 ■ p 2 2p 1 ■ p 2 \ 

2 \ (ki ■ k 2 )(ki - pi) (h ■ k 2 )(h ■ p 2 ) (h ■ p^ih ■ p [ 2 JJ 

As in sect. 3, we fix 9 = k2Pl , and integrate over k\ (k 2 < k\ < pi) to obtain 
the 1PI cross section. 

At a large angle (0/9 piP2 = 0(1)), the condition 9k 1 k 2 ^ @k 2 p2 — @k iP2 or 
9k ipi 9 P1P2 — @k iP2 is n °t necessarily satisfied, and we cannot make use of 
the approximation (20). We therefore have to analyse each term in (36) more 
carefully. 

The /ci-integration of each term can be performed with the use of the 
decomposition similar to that in subsect.2.1. The azimuthal integration of 
each term gives 

k 2 ■ p x k\ 1 [ 1 1 1 



(h • k 2 )(k 1 -pi) 2 { (1 - cos6» fclfe2 )(l - cos6» fclPl ) l-cos0 felPl l-cos6> fclfe2 

+(pi <-»■ h) 



d(cos9 klP1 ) 

"©(^fcaPi - d k lPl ) + (pi <-> k 2 ), 



h(l- cos9 klPl ) 



k 2 ■ p 2 k\ 



+ 



(fci • k 2 ){ki -p 2 ) 2 { (1 - cos6» fclfc2 )(l - cos6» felP2 ) l-cos6> felP2 l-cos0 fclfe2 

+ (p 2 k 2 ) 
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fci(l - cos0 klP2 ) 

Pl-P2fc? _lf 1 | 1 1 1 

(ki ■ pi)(ki ■ p 2 ) 2 \ (1 - cos0 fclP1 )(l - cos9 klP2 ) l-cos9 klPl l-cos6 klP2 j 

+{pi <-> p 2 ) J 1 71 — ^~\tt — v + r~ r" - i + (pi <-> P2) 

We thus obtain 

C^g s d 3 A:i / 2k 2 ■ pi 2k 2 ■ p 2 2gi ■ p 2 | 

(2tt) 2 A;i \(A;i -fc 2 )(A;i-pi) (h ■ k 2 )(h ■ p 2 ) (h ■ p l ){k 1 ■ p 2 )\ 

_l_ /* ^pJ_ [@(^ 2p2 _ 6> fcip2 ) - &(9 PlP2 - 9 kiP2 )\ 

J 1 — COS t7fc 1 p 2 

+ / d(cos y [9(0 - fclfc2 ) + 6(0 fc2P2 - felfc2 )]| (38) 

J l-COStf fcl fc 2 J 

Note that only the last integral (over cos^^) on the rhs of (38) gives 
logarithmic divergence. Let us isolate the divergent part in the form of the 
exact AO 

J l-cos0 felfc2 
The rest gives its correction: 

r cosd d(cos6 klPl ) _ r^e k2P2 d(cos0 fclP2 ) /" cose d(cos0 fclfc2 ) 

Jcos9 P1P2 1 — COS 9 klPl Jcos9 P1P2 1 COS 9 klP2 Jcos6k 2P2 1 — COs9 kl k 2 

In ^ cos 9 PlP2 \yl^ cos 9 P1P2 ^ 1 cos 9 k2P2 

1 — cos 1 — cos 9 k2P2 1 — cos 

= ^ln 1 " 005 ^ 2 . (40) 
1 - cos 9 k2P2 
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Particularly, in the centre-of-mass frame of the gg-pair (9 PlP2 = ir, 9k 2P2 = 
n 0). (40) is 

- 2 In — n = -2 In = 2 In [ cos 2 -J . (41) 

1 - cos 9 k . 2P2 1 + cos 9 \ 2 J v ' 

Note that (41) vanishes as 9 — > 0. 

Now let us generalise the above result to all orders. At the next-to-leading 
order, we can make use of (34): We identify the direction of the registered 
particle with that of hi, and multiply the multiplicity of the /^-subjet to 
the A; 2 -gluon emission cross section. The opening angle of the /c2-subjet is 
identified with 9k 2 ki, owing to the AO. 

We calculate the correction to the AO with the use of solution of the fixed 
coupling integral equation (29). In the result, we can replace the solution 
with running coupling (32) for M g , and the value of 7 (33) for the anomalous 
dimension. Indeed, we can calculate the correction using running coupling 
from the start without any difficulties. Because the difference appears only 
at the next-to-next order, we use the fixed coupling for the sake of notational 
simplicity. 

When we discussed the resummation in the previous section, all emission 
angles are assumed small. If they are not necessarily small, we have to be 
careful of the choice of the variables. The emission angle, if it is large, is 
not invariant under the boost in the direction of the parent parton. In fact, 
the resummation is organised in terms of the transverse momentum relative 
to the parent parton. The multiplicity of the /c 2 -subjet (with the fe-gluon 
emitted from the fci-gluon) is given by 

M g = (k 2 2 sin 2 9 klk2 y (42) 

where k 2 sin# felfe2 is the transverse momentum relative to k±. In (42) we have 
put unity for the normalisation constant C in (29) for notational simplicity. 
With the AO part of the integration kernel (39), we obtain 

dcOS^fa Mg{k 2 shl 2 dkik2)e{d _ dkik2) 

-- - (2/c|(l - cos #)) 7 . (43) 
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This result is not written in terms of the transverse momentum k2± = 
k2 sin 9, and not invariant under the boost in the jet direction. 

Let us now add lhs of (40) to the integration kernel on the lhs of (43). 
In the correction, the integration kernel over the angular variable does not 
become singular in the integration range. It implies that the angular variable 
in the argument of the multiplicity function may be set equal to 9 (= 0(1)). 
(We may put 9k 2P2 or # fclfc2 in place of 9. In the absence of the singularity in 
the kernel in the integration range, the difference in the result is at higher 
orders.) 

- r " ^2!%j) M ,(kl S ,n> S) - r" hm '" {C " S \" ) M s (kl sin* S) 

Jcos0p lP2 1 — COS Ui si p 1 J cos Sp lP2 1 — COS t/fc 1 p2 

Jcos8 k2P2 1 — COS f/fc 1 fe 2 

= 2 In ^cos 2 ^ {k\ sin 2 9)"' + higher orders . (44) 

Adding the correction (44) to (43), we obtain 

^ (±k 2 2 sin 2 ^ 7 + 2 In (ms 2 ^ (k 2 2 sin 2 9f 
2 

= -{k\ sin 2 9) 1 + higher orders . (45) 

7 

The result is now boost invariant. 

The total multiplicity is given by integrating (45) over the energy variables 
k 2 , ki, and the angle variable 9. The 9 integration of the AO part, the rhs of 
(43), with the integration kernel 1/(1 — cos 9) gives (2/7 2 )(2/c|) 7 . The similar 
integration of the correction, the rhs of (44) over 9, on the other hand, gives 



1 d(cos6») 2 2 . 2 

- / — i 7r21n- -(/^sin 2 ^) 

J-i 1 - cos 9 1 + cos6» v 2 y 



= -2C(2) (A; 2 ) 7 + higher orders , (46) 

where 

. , v f 1 d^ , 1 7T 2 , 

C(2 = / . ln - = 7T- ^ 

JO 1 — z z o 
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The correction, therefore, reduces the total multiplicity obtained with the 
AO by a factor (1 — 7 2 7r 2 /6). This result agrees with the dipole correction 
calculated by Dokshitzer et a£.||0 using the planar gauge. 



6 Conclusion 

In this article, we have examined the next-to-leading order corrections to the 
AO in the angular distribution of particles in a jet. At the leading order, 
the angle of the registered particle is identified with the direction of the 
subjet to which it belongs: i.e. that of the primary gluon, which is directly 
emitted from the initial energetic parton. At a small angle (i.e. near the jet 
direction), we found that there is no next-to-leading corrections to the AO. 

In the particle flow generated by the coherent Monte Carlo simulation, 
the large-angle emission of a primary gluon contributes to the forward mul- 
tiplicity, while such contribution is absent in the exact cross section. On the 
other hand, the contribution from the small-angle emission of the primary 
gluon is somewhat reduced, and the angular distribution obtained in this 
simulation is correct at the next-to-leading order. 

At a large angle, the correction to the AO reduces the multiplicity by 
a factor of (1 + 7 In cos 2 (#/2)). The correction restores the boost invari- 
ance of the angular distribution, which is absent in the AO formalism. In 
another word, the dipole correction appears only because we took the an- 
gle variables to formulate the coherence. There is no such correction at the 
next-to-leading order, if we use Lorentz-invariants (in place of angles) in the 



formulation [ 15 1 [ 13 



The value of the anomalous dimension 7 is around 0.20 for a s ~ 0.12 
(around the LEP energies). The correction at 6 = n/3, for example, is about 
5.7%. It is comparable to size of the uncertainties from the yet uncalculated 
higher order corrections. 
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FIGURE CAPTIONS 



Fig.l The lowest order diagrams for one-gluon emission cross section. 

Fig. 2 The diag ram for the soft-gluon emission (at a momentum A^) from 
a narrow jet of q + g. We assume &2 • Pi <C &i ■ p±. 

Fig. 3 Four angular regions in the AO approximation. The initial energetic 
qq is produced at the rectangle. There is no soft-gluon emission in the region 
VI. 
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